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Abstract 

The lowest order resonant bifurcations of a periodic orbit of a Hamiltonian sys- 
tem with two degrees of freedom have frequency ratio 1 : 1 (saddle-centre) and 1 : 2 
(period-doubling). The twist, which is the derivative of the rotation number with 
respect to the action, is studied near these bifurcations. When the twist vanishes 
the nondegeneracy condition of the (isoenergetic) KAM theorem is not satisfied, with 
interesting consequences for the dynamics. We show that near the saddle-centre bi- 
furcation the twist always vanishes. At this bifurcation a "twistless" torus is created, 
when the resonance is passed. The twistless torus replaces the colliding periodic orbits 
in phase space. We explicitly derive the position of the twistless torus depending on 
the resonance parameter, and show that the shape of this curve is universal. For the 
period doubling bifurcation the situation is different. Here we show that the twist 
does not vanish in a neighborhood of the bifurcation. 

Keywords: Twist Maps; Hamiltonian Systems; Saddle-Centre Bifurcation; Period- 
doubling Bifurcation; KAM; Normal Forms; Elliptic Integrals 

1 Introduction 

The dynamics near a periodic orbit of a Hamiltonian system can be studied in terms of a 
local Poincare section transversal to the orbit. In two degrees of freedom the first return 
map restricted to the surface of constant energy is an area preserving map with a fixed 
(or periodic) point. If the multipliers fa, i = 1,2 of the fixed point have modulus 1 (but 
are not equal to ±1) the fixed point is called elliptic. Then [i\ = /ZJ G C and we can 
write [i = exp(27riw) with the rotation number < u < 1/2 of the periodic orbit. If the 
periodic orbit is elliptic and the rotation number is irrational the map can (formally) be 
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transformed to Birkhoff normal form which in action-angle variables (<p, I) reads 

(^7)^(^ + 27^(7),/). (1) 

The action 7 is like the radial coordinate in polar coordinates, hence the map in normal 
form maps circles to circles by rotating them Sl(7) times. The rotation number (or winding 
number) £1(1) near the periodic orbit can be expanded as 

Q(I) =lo + t I+^t 1 I 2 + ... 

The twist (or torsion) r(7) is the derivative of the rotation number with respect to the 
action, 

t(I) = ^(I)=t + t 1 I+... 

When the rotation number is a strictly monotone function of the action in some interval the 
map (^Q) restricted to the corresponding annulus is called a monotone twist map. Moser's 
KAM theorem states that the invariant torus 7 = Iq of Q persist under perturbation 
when its frequency £l(Io) is diophantine and its twist r(7o) does not vanish. Arnold's 
KAM theorem [Q is the same statement for flows where the nonvanishing of the twist 
corresponds to the isoeneregetic nondegeneracy conditions. A well known corollary of the 
KAM theorem is the stability of an elliptic fixed point in two degrees of freedom when 
O(0) = w/ 1/3, 1/4 and the twist at the origin is non-vanishing, r(0) = tq ^ 0. 

When the twist vanishes the perturbed dynamics can be more complicated. The sta- 
bility of an elliptic point can be lost when its twist vanishes, see p] for an example of an 
unstable elliptic point with uj = 1/5. The effects of vanishing twist away from the origin 
was first described by Howard [Jj, and the resulting effects have been observed in many 
examples ^JJE3- The probably most spectacular effect is the appearance of so-called 
meandering curves H3 El • The properties of non- twist maps also show interesting be- 
haviour under renormalisation [2] and recently it has been shown [3| that an extension 
of the KAM theorem can also be proved in this context. In jSJ it was finally shown 
that the vanishing of twist at the fixed point generically occurs in a one parameter family 
when the rotation number of a fixed point passes through the interval [1/4, 1/3]. When 
the twist vanishes at the fixed point a twistless torus is created in a twistless bifurcation 
[3] . After creation the twistless torus passes through resonances and in this way non- twist 
maps generically appear in one parameter families of area preserving maps. The truncated 
resonant Birkhoff normal form shows that this twistless torus eventually collides with a 
saddle-centre bifurcation that gives rise to the period 3 orbits that collide with the fixed 
point when u = 1/3 5 . Such a connection between resonance and vanishing twist can 
also be found in 4 dimensional symplectic maps [3]. 

The techniques of j^j can also be applied to the higher order resonances, in particular 
for u = 1/4. In this paper we study the two remaining generic bifurcations at even 
stronger resonance uj = 0, 1/2. By definition the corresponding fixed point is not elliptic. 
The two cases will be denoted as the 1 : 1 and 1 : 2 resonance, or as the saddle-centre and 
period doubling bifurcation, respectively. The main result is that near the saddle-centre 
bifurcation the twist always vanishes, while it does not vanish near the period doubling 
bifurcation. 

The method is based on the analysis of the resonant normal form, in which the Poincare 
map is approximated by the time 1 map of a one degree of freedom system, see e.g. 8 . This 
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Figure 1: Lines of constant energy on the phase space (u,v) for the saddle-centre bifurca- 
tion (1 : 1 resonance, to = 0). left: e = —0.1, middle: e = 0, right: e = 0.1. 

normal form is an approximation, that is local near the bifurcation in parameter space and 
in phase space. At first we will completely ignore this, and just analyse the normal forms 
in the following two sections. In section 0] we will address the problem of non-locality in 
phase space, and also comment on the effect of higher order perturbations on the twistless 
tori. We will talk of invariant tori even though the invariant curves H(u, v) = h of the 
normal form may not be compact. This will also be justified in Sec. El Finally we treat a 
saddle-centre bifurcation in the Henon map as an example. 

2 Saddle-Centre Bifurcation 

The normal form of a Hamiltonian system with two degrees of freedom near the 1 : 1 
resonance has the form 

iP 1 up 

H( u ,v) = — + — + eu. (2) 

The coefficient of it 3 has been scaled so that it equals 1/3. The variables u and v are 
canonically conjugate variables on a local transversal Poincare section and e is a parameter, 
typically the energy of the original system. The Poincare map is given by the time 1 map 
of the flow of H(u, v). If the period T is large, the time 1 map advances little. The rotation 
number of the full system is the period 1 divided by the period of the reduced one degree 
of freedom flow, f2 = 1/T. The more familiar Q = 2tt/T is obtained when the time 2tt 
map is taken instead of the time 1 map, but the time 1 map is more natural at least for 
the example of the Henon map we are going to discuss. Since f2 is determined by T, we 
now study in detail the period T of the one degree of freedom system given by H. 

The critical points and critical values of the energy map H : R 2 -► R and their 
dependence on e give the main structure to the bifurcation. Instead of a one degree of 
freedom system H(u, v; e) depending on the parameter e one may consider H(u, v, e, 9) 
as a Hamiltonian in R 3 x S , with action e and conjugate angle 9. The set of critical 
values of the energy-momentum map (H , e) : R 3 x S 1 — > R 2 is called the bifurcation 
diagram. A simple way to compute it is to find the critical values of the energy map 
of H(u, v) and consider their parameter dependence on e. The Hamiltonian has critical 
points (u, v) = (±-\/ = e, 0) and corresponding critical values h = =F2(— e) 3//2 /3. They exist 
when e < and the upper sign corresponds to a local minimum of H, while the lower sign 
gives a saddle. The corresponding phase portraits are shown in Fig. ^ 

The dynamics is given by Hamiltons equation u = v and eliminating v using the 
Hamiltonian gives a first order equation for u. After separation of variables the period of 
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motion with energy h is given by the elliptic integral 



' P , dU (3) 

'2h - |n 3 - 2eu 

where the integration path is encircling the interval on the real axis where the argument 
of the square root is positive. If there are two positive intervals either one can be taken, 
the result is the same. By scaling u = z(ae) l l 2 , where a = sign(e) and introducing the 
one essential parameter 

3k (4) 



' 2(ae) 3 / 2 ' 

the period T is an elliptic integral on the curve 

£ : w 2 = P 3 (z) = 2 7 - z 3 -3az. (5) 

The case e = has to be excluded in this scaling, but it is simple to treat it separately. 
We are mostly interested in the case where a = 1. The essential integral now reads 

S{ 1 ) = j^dz (6) 

and it is related to the period by 



r(/l ' £) = (^ 5(7) - (7) 

The polynomial P3 has one or three real roots. The collision of two real roots corresponds 
to the unstable equilibrium and its separatrix. It occurs when the discriminant 

A = -108(o- + 7 2 ) 

vanishes. A = is only possible for a = — 1 and hence the critical parameters for which a 
double root occurs are given by 7 = ±1, hence 

9h 2 = -Ae 3 , (8) 

which has a cusp at the origin. At the origin h = e = all three roots collide in the 
saddle-centre bifurcation. The discriminant (jSJ) is shown in Fig. [21 In the case of the 
saddle-centre bifurcation the bifurcation diagram is given by the discriminant of P3. The 
bifurcation diagram divides the parameter plane (e, h) into two regions: one with 3 real 
roots and one with 1 real and two complex roots. The latter has positive e everywhere, 
while the former is the wedge shaped region contained in the negative half-plane. For 
e < the phase portraits contain a pair of stable/unstable fixed points. The critical value 
of the energy of the stable fixed point is a local minimum given by the lower branch of the 
bifurcation diagram, while that of the unstable fixed point is a saddle given by the branch 
with positive h. For this range of energies there are 3 real roots. It will turn out that 
the upper branch of the bifurcation diagram corresponding to the unstable fixed point 
is crucial for the existence of vanishing twist. For e > the phase portrait is without 
fixed points. Even though the topology is trivial in this case we will now show that the 
rotation number has a maximum on a certain invariant torus containing points near the 
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Figure 2: Schematic sketch of the bifurcation diagram for the saddle-centre bifurcation. 
Graphs of P3 are shown together with a horizontal line indicating the value of h. The bold 
lines are the critical values of the saddle. 




Figure 3: Lines of constant period T = 1/Q equidistant with AT = 0.3 on the parameter 
plane (e, h) for the saddle-centre bifurcation (1 : 1 resonance, lu = 0). The vanishing twist 
is indicated by a curve of vertical tangents dil/dh = 0. 
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origin in the phase space. At this maximum of the rotation number the twist vanishes. 
The vanishing twist occurs at the vertical tangents of the contours of the rotation number 
shown in Fig. EI The fact that the invariant curves are all unbounded for e > will be dealt 
with in section 0J For now observe that the integral © is finite, even though the invariant 
curves are unbounded in v. The main feature of the level lines of the rotation number 
as shown in Fig. |3] is that it diverges when the unstable periodic orbit is approached. 
This occurs for the positive critical value of h when e < 0. Everywhere else the rotation 
number is a well defined, smooth and bounded function of h and e. Accordingly the level 
lines "hug" the curve of critical values that correspond to the unstable orbit. Already 
from this property one can deduce the existence of a curve with vanishing twist using 
topological arguments. Here we proceed along the analytical route, because it will give 
us more detailed information. Note that the curve of critical values with negative h, see 
Fig. does not appear in Fig. |21 One reason for this is that we chose to plot the rotation 
number for the (non-compact) invariant tori with motion between — oo and the smallest 
real zero of P3. These invariant tori do not contain critical points of the energy map, 
even though the corresponding energy might be a critical value. Accordingly the rotation 
number is smooth across this line of critical values. The critical point corresponding to 
the critical values is the stable fixed point at the local minimum of H. But even if we 
would plot the rotation number of the bounded invariant tori near the local minimum of 
the potential the picture is unchanged. The reason is that for a cubic elliptic curve the 
integrals of first kind over either one of the two real intervals (if they exist) are equal. 
Since f2 = 1/T and 2irT = dJ/dh the derivative of the rotation function is 

dn _ m/dh dT/dh 
~dJ ~ dJ/dh ~ ~ 2nT 3 

Hence the twist vanishes when 

&T _ (3/2) 3 / 2 OS 

dh ~ £7/4 ~ 

and this is only possible for finite e when 

dS / 1 

— = - j)—dz = 0. 

07 J ttr 

This complete elliptic integral can be written as a linear combination of Legendre's stan- 
dard integrals. In this way a condition for the vanishing of the twist is now obtained. The 
relevant case for this purpose is that of one real root, for which the phase portrait has no 
fixed point. The integrand w is positive for u S (—00,20)5 where zq is the single root of 
Ps(z). Let the factorized polynomial be given by 

P 3 (z) = -(z-z )P 2 (z), J p 2 (z) = ( z - Cl ) 2 + C 2 2 , (9) 

so that the complex roots are Ci ± K2- Denote the distance between the real and complex 
roots by r, hence r 2 = /^(^o), so that the discriminant is given by A = — 4r £|. Legendre's 
standard integral of the first kind K(k) has differential 

ui = - 1 dz (10) 
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up to a constant factor, where the modulus k is given by 



2 V r 7 2 V v'TT 



Q 

and in the last equality the parameter a = £2/(^0 ~~ Ci) nas been introduced. In the second 
equality in addition sign(zo — Ci) = sign(7) is used, which is true because zq = Q\ in Q 
together with the vanishing of the quadratic coeffcient in Q implies zo = Ci = 0) an d 
therefoe the polynomial has no constant term and h = 7 = 0. 

A non-standard form of the differential of Legendre's standard integral of second kind 
E(k) is 

P 2 (z)dz 
(z - [z ± r)) z 

up to the same constant factor as in (|1U|) . The differential dz/w 3 we are interested in is of 
the second kind, and can therefore be written as a linear combination of u\ and u>2 with 
constant coefficients, up to a total differential: 

^ = Auji + Buj 2 + dF 
ur 

where F = Q2{z)/{{z — (zo ± r))io). Together with the undetermined coefficients of the 
quadratic polynomial Q2 this gives a system of 5 linear equations for the 5 unknown 
coefficients. Solving these equations gives 

A = ^((zo - (i)r - r 2 + 4(1), B = —(r 2 — 4C|) 

Since the quadratic coefficient of P3 is zero, the roots of P3 add up to zero. Therefore 
the real parts satisfy zq + 2£i = 0, hence 

2/c 2 = I-3Ci/r and r 2 = 9(f + Cf . 

With these equations the coefficients A and -B can be expressed in terms of k alone, up to 
the factor (r 2 A) _1 . The condition of vanishing twist, dS/d'y = 0, finally reads 

(8k 4 - 9k 2 + l)K(k) = {16k A - 16k 2 + l)E(k) , 

where K(k) and E(k) stand for elliptic integrals of the first and the second kind, respec- 
tively. 

For k = we have equality since both elliptic integrals equal tt/2. The first derivatives 
of either side vanishes, but the second derivatives are — 35-7r/4 and — 657r/4, respectively, 
so that the left hand side is larger for small k. For k = 1 the prefactor of K vanishes, 
while that of E gives 1 and E(l) = 1. Hence for k — > 1 the right hand side dominates. 
This proves that there exists a solution of this equation for k S (0,1). Numerically we 
find k 2 S3 0.7097215. In order to calculate the corresponding 70 we observe that a as 
introduced in (fTT|) is related to 7 by 

C2 C 2 1 r 2 +a 
a -z -Ci" 3Ci-V3~r 2 -a' [L2) 

where 



r = (7 + v^+W 73 • (13) 
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Figure 4: Bifurcation diagram of the position in phase space u versus the bifurcation 
parameter s. For e < the location of the fixed points at ±V — e are shown, while for 
e > the maximal u of the twistless torus at ^o(7o)\/£ is shown. 



Using and ko the corresponding value of a is ao ~ 2.164255, and from ao using (|12|) 
we find 70 ~ 0.9152203. Therefore the curve of vanishing twist in the parameter plane 
occurs for positive e when 

3h = 2 7o £ 3/2 . (14) 

Since 70 < 1 the curve of vanishing twist is bent downward as compared to the bifurcation 
curve (JBJ) for e < and h > 0. See Fig. |3] for a graph of this curve together with the 
numerically computed lines of constant rotation number. The lines of constant rotation 
number have vertical slope at their intersection with the critical curve, as must be the 
case. 

The scaling of u reduces the number of parameters to one. The essential parameter 7, 
in its dependence on h and e, organizes the bifurcation. It allows to compute explicitly all 
the important characteristics of the twistless torus. Combining (jllj) and (|12j) shows that 
A; is a function of 7. The curves in the parameter plane that have the same value of 7 are 
given by (|4jl. The most prominent ones are the curve of critical values 7 = ±1, as shown 
in Fig. |21 and 7 = 70, the curve of twistless tori, see Fig. |31 They all have the same shape 
of a semicubical parabola, except when 7 = 0, hence h = 0, or 7 = ±00, hence e = with 
±h > 0. 

k 2 7 a curve 



+1 0^ 3/i = -2(-e) 3 / 2 

i(2 + V^) +00 l/x/3 e = 0,h>0 

0.709721497 0.91522 2.164255 3h = 7o(-e) 3/2 
\ 00 e = 0,h > 

i(2-V3) -00 l/y/3 e = 0,h<0 

-10 3/i = 2(-e) 3 / 2 

Approaching the bifurcation point along the curve 1(3} gives ^(7) in the limit. The function 
k(h, e) is therefore not continuous at the origin. Moving on curves (JJJ) in the parameter 
plane for any 7 / 1 the change in the period T(h, e) given by (jHJ) is elementary. From (J2J) it 
follows that T is proportional to |e| -1 / 4 , and the constant of proportionality is determined 
from (J7J) and ©• The divergence of the period upon approaching the bifurcation is 
therefore not caused by the elliptic integral, but merely by the algebraic dependence 



S 



|e| 1 / 4 . In this way the rotation number of the stable periodic orbit is given by 
Q( 7 = -l, e ) = -^(-e) 1 / 4 « 0.225079 \e\ 1 ^ . 

V 27T 

The integral S(—l) contained in this expression can be easily calculated using residue 
calculus because for 7 = — 1 the curve has a double root, P%{z) = — {z + 2)(z — l) 2 . For 
general values of 7 and in particular for the twistless torus the elliptic integral S(j) © 
needs to be calculated. The single real root is given by 

z = r - p hence z (7o) ~ 0.5535942 . (15) 

Finally 5(7) is obtained as 
4K(k) 



'3z \/l + a 2 /2 

In particular when 7 = 70 the rotation number of the twistless torus is 

The constant of proportionality is close to v3/47r. Using the above value of zq the position 
of the rightmost point of the twistless torus in phase space is located at «o = z^^fe. This 
means that for e = £q this point is on the same side as the stable periodic orbit was for 
£ = —£o before the bifurcation, but by a factor of l/zo(7o) ~ 2 closer to the origin. See 
Fig. 0] for an illustration in the form of a standard bifurcation diagram showing position 
in phase space u versus bifurcation parameter e. 



3 Period-Doubling Bifurcation 

The normal form of a Hamiltonian system with two degrees of freedom near a periodic 
orbit in 1 : 2 resonance is 

H(u,v) = y + D^ + eu 2 , (17) 

where D = ±1. As in the case of the saddle-centre bifurcation the variables u and v are 
canonically conjugate variables on a local transversal Poincare section and e is a param- 
eter, typically corresponding to the energy of the original system. Contourplots of this 
Hamiltonian show the intersection of invariant tori with the Poincare section transversal 
to the bifurcating orbit. For D = 1 they are shown in Fig. for D = — 1 in Fig. El for 
e < 0, e = 0, and e > 0, respectively. 

The critical points and critical values of the energy map H : R 2 R and their 
dependence on the parameter e describe the structure of the bifurcation. The Hamiltonian 
has critical points at (0, 0) and at (u, v) = (±V — 2eL>, 0) with critical values h = and 
h = —D£ 2 . The origin is a local minimum for e > 0, a saddle otherwise. The second 
critical point exists when sD < 0, and is a local minimum for e < 0, a saddle otherwise. 
The set of critical values, see Fig. [71 therefore is the union of the line h = with the half 
of the parabola h = —D£ 2 for which De < 0. For D > the only unstable branch in the 
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Figure 6: Lines of the constant energy on the phase space (u, v) for the period-doubling 
bifurcation (D = —1), for e = —0.1, e = 0, e = 0.1 



bifurcation diagram is h = for e < 0. It divides the two regions of real motion. A third 
region {h < —£ 2 } U {e > 0, —e 2 < h < 0} is not in the image of H , hence there is no real 
motion corresponding to (h, e) from this region. 

For D < the line of critical values /i = again has a saddle as critical point when 
e < 0. In addition the half parabola h = e 2 , e > also corresponds to a saddle of H that 
is not at the origin. When D < the whole plane (h, e) is in the image of H considering 
all e. The critical values divide the plane into three regions with 0, 2, and 4 real roots. 

The dynamics of the reduced one degree of freedom system is given by ii = v and 
eliminating v using the Hamiltonian gives a first order equation for u. Separation of 
variables then gives the period of motion in the reduced one degree of freedom system as 
the elliptic integral 

du 

^2h - \Du A - leu 2 

The number of parameters could be reduced by introducing the ratio h/e 2 , but for clarity 
we do not introduce this scaling. The period T is an elliptic integral on the curve 

£ : w 2 = P 4 (z) = P 2 (z 2 ) = 2h- ^Dz 4 - 2ez 2 . (18) 

The discriminant of P4 is 

A = -256Dh(Dh + e 2 ) 2 . 

It vanishes for h = and Dh = —e 2 , corresponding to the critical values h = and 
h = -De 2 already found above. As usual the set of critical values is contained in the 
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D=1 



D=-1 




Figure 7: Schematic sketch of the bifurcation diagram for the two cases of the period 
doubling bifurcation. Graphs of P4 are shown together with a horizontal line indicating 
the value of h. The bold line denote critical values of unstable orbits. The dashed line is 
not critical; it indicates a vanishing of the discrimiant. 



discriminant, however, the discriminant vanishes on a larger set. Namely the branch of 
the parabola with De > is not part of the critical values. The three regions already 
found correspond to regions with 4, 2, and real roots of P4, respectively, see Fig. d The 
part of the discriminant that is not part of the critical values is dashed. The polynomial 
can be factored as 



P 4 (z) 



1 



1 



D(^-z 2 )(Z + -z 2 ), -DZ± = -e±Ve 2 + Dh. 



Comparing coefficients gives = —4hD and +£+ = — 4eD. The factorization of P4 
has real factors, i.e. £± £ R, with one exception. It occurs when the quadratic equation 
P^iO = has complex roots. The position of the roots in the regions of the parameter 
plane is as follows. The two real cases for D = 1 are 



+2 : h > 

+4 : e < 0, -e 2 < h< 



£- < < 
< £- < 



2[-ve 



In the first column a label is given containing the number of real roots. All 4 cases appear 
when D = —1 for 



-2 
-4 
-Oi 
-0c 



/»< 

e > 0,0 < h < e 2 
e < 0,0 < h < e 2 
h>e 2 



e+ < < e_ 
< e+ < e_ 
e+ < e- < 

£± complex 



2[-oo,-ve-] 

2[-oo,-7^] 
[—00, 00] 
\— 00, 00 1 



The last column gives the interval of real motion along which the period integral is taken. 
The region without real roots contains two parts separated by the branch of the parabola 
h = —De 2 with De > 0. On this branch there occurs a collision of complex roots at 
z 2 = —2eD, and they move from the imaginary axis into the complex plane. 

The intervals of integration and multiplication factors as given in the last column of 
the previous table can be read off from the phase diagrams, see Fig. El In the case 
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Figure 8: Lines of constant period T = 1/Q equidistant with AT = 0.3 on the parameter 
plane (e, h) for the period-doubling bifurcation (1 : 2 resonance, u = 1/2) with D = 1 



D = — 1 with 4 real roots there are different orbits for the same (/i, e). One is the compact 
orbit near the stable fixed point with extent 2[— y^+, v/£+]- The interval given above is 
for the non-compact orbit, which has half of this period. 

In order to calculate the derivative of T with respect to h the integral is first written 
in standard form and then differentiated. This more traditional approach (as compared to 
the previous section) is preferable in this case because the roots of the quartic are easily 
written down. 

Denote the ratio of the roots £± by 



_ £_ _ -e - Ve 2 + Dh 
r ~ £+ ~ -e + Ve 2 + Dh ' 
Then the period T(h, e) in the 6 cases is given by 

+2,-2: T= / 8(l-2fc 2 ) \ 1/2 ^ fc 2 = -J-, (19) 



+ I. -()/ : T - I 8(fc2 2) ) ' K(k), k 2 = l-r, (20) 



i: r-[^-^l) 1/2 K(k), fc 2 = i, (21) 



The level lines of the period T (and hence the rotation number 0, = 1/T) are shown in 
Fig. [HI for D = 1 and in Fig. El for D = — 1. These numerically computed pictures show 
that there are no vertical tangents, hence the twist does not vanish. This is now proved 
by differentiating the period in Legendre normal form. 
The twist vanishes when 

_ dT _ dTdk 
dh dk dh 
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Figure 9: Lines of constant period T = equidistant with AT = 0.3 on the parameter 
plane (e, h) for the period-doubling bifurcation (1:2 resonance, u> = 1/2) with D = — 1 



The last factor dk/dh does not vanish for k 6 (0,1). In the cases +2, —2, and — 0c it 
seems to vanish when k = 1/2. But this implies that e = and the singularity cancels. 
Therefore it is enough to consider dT/dk = 0. After removing common non- vanishing 
factors the conditions for vanishing twist are Vi(k) = 0, i = 1,2,3 where 

+2,-2,-0c: Vi{k) =(1 - 2k 2 )E(k) - (I - k 2 )K(k), (23) 

+4, -Oi : V 2 (k) = (i " y) E ( k ) ~ (1 - (24) 

-4: F 3 (fe) =(l + k 2 )E(k)-(l-k 2 )K(k). (25) 

The three equalities Vi(k) = 0, i = 1,2,3 are never satisfied on the range k € (0,1). 
Obviously V$(0) = 0, while Vi(l) = —1,1/2,2 for % = 1,2,3, respectively. So we need to 
show that Vi is negative and V2, V3 are positive for A; G (0,1). Differentiating Vi and V2 
gives the simple results 

Vi{k) = - h(K(k) - E(k)) (26) 
V^'(jfe) = - 3kE(k) (27) 

In the last case V^(k) is obviously non-positive, so that the twist in the case —4 is a 
monotone funtion rising from 0, hence it is nonzero. Similarly also V^ik) is non-positive, 
which follows from the well known inequality K(k) > E(k). The first function V\{k) is 
negative on (0, 1), but not monotone. Rewriting it as 

Vi(k) = (1 - k 2 ){E(k) - K{k)) - k 2 E(k) < 

the inequality is clear because both terms are negative for k £ (0, 1). Therefore the twist 
never vanishes in a neighborhood of the period doubling bifurcation. 

The relation between the cases D = ±1 is interesting. The main observation is that 
changing the sign of D and e inverts the overall sign of the potential V(x; D, e) = Dx 4 /4 + 
2ex 2 . Therefore changing the signs of D, e, and h leaves the roots of P4 invariant. By this 
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mapping the regions in parameter space with the same numbers of real roots are mapped 
into each other. The integration path does change in a less trivial way. The integration 
needs to be taken over the positive intervals of P4 on the real axis. Changing the sings 
of D, e, and h does change the sign of P4. As a result the periods for D = 1 are the a 
cycles of the elliptic curve, while those for D = — 1 are the (3 cycles. Hence the period for 
the case —4 can be obtained from that of +4 by replacing K{k) by K'(k) = K(\/l — k 2 ). 
Similarly, the cases +2 and —2 are mapped into each other. 

The obvious symmetry in Figures |H] and El with respect to changing the sign of h is 
related to the fact that changing the sign of D and h leaves the ratio r, and therefore also 
the corresponding modulus k 2 , invariant. This means that the level lines in region +2 and 
—2 can be obtained from each other by reflection on the h = axis. In a similar way the 
regions +4 and — Oi have the same rotation number. 

The period doubling could have been treated in a scaled version with only one essential 
parameter 5 = h/e 2 . However, the presentation seems more transparent in the unsealed 
version. Similar to the case of the saddle-centre bifurcation the modulus k 2 of the elliptic 
integral for the period is constant on the parabolas h = 5e 2 . Again dependence on the 
parameters on these curves is simply algebraic, as before through |e| 1//4 . The value of 
the modulus is not defined at the origin, but depends on the parabola on which it is 
approached. But in any case, there are no twistless tori near the origin in this bifurcation. 

4 Universality 

A major problem in our approach seems to be that the invariant twistless tori in the 
normal form are not compact. The normal form is obtained from an expansion near 
the bifurcation point, and is therefore local in phase space and local in the parameter. 
How can the rotation number of a non-local invariant torus be determined from this local 
normal form? To answer this question higher order terms need to be considered in the 
Hamiltonian. The Poincare map near the bifurcation can be described by the time 1 map 
of the Hamiltonian 

H(u, v, t; e) = H(u, v; e) + G(u; e) + R(u, v, t; e) . 

The remainder terms R containing the periodic time dependence can be made arbitrarily 
small, but they cannot in general be removed while retaining a non-zero radius of conver- 
gence of the normal form. The first term H is the normal form analysed in the previous 
chapters. We will concentrate on the saddle-centre case (j2J), but similar remarks apply 
to (|17|) , The invariant tori near the origin for e > of H are not compact. The higher 
order terms in G can compactify them. The results about vanishing twist can be applied 
when G does compactify these curves. However, the precise form of G does not matter. 
Under the compactness assumption KAM theory can be applied to H + G, where R is the 
perturbation. Many of the invariant tori of H + G will persist. In particular a twistless 
invariant torus of H + G will persist if it is sufficiently irrational. The curve (|14|) in the 
parameter plane therefore does not have invariant twistless tori in its perimage for every 
point, instead just for a cantorset of points. This is well understood. The main issue 
in this section is to understand the effect of adding G to H. For small e the essential 
contribution to the diverging period comes from the dynamics near the origin, while the 
dynamics on the invariant torus away from the origin has finite period. This is why the 
local normal form can give a statement about the dynamics on a non-local invariant torus 
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near the bifurcation point. In particular we will now show that the curve of vanishing 
twist that was found to be emanating from the cusp of the saddle-centre bifurcation has a 
universal shape sufficiently close to the cusp singularity. In particular this means that the 
constant 70 that determines the shape of the curve of twistless tori in relation to the curve 
of critical values of the unstable orbits has the universal value 70 ~ 0.91522. Quantities 
derived from 70, like zq and the coefficients in (|16|) . are accordingly also universal. Our 
calculation will show that the value of 70 is not influenced by the higher order terms G 
in the Hamiltonian. Moreover, the following argument also shows that integrating the 
non-compact invariant torus up to v = ±00 does not introduce and additional error. 
Let the high-order truncated normal form Hamiltonian be 

H(u, v; e) = \-v 2 + ^n 3 + eu + G(u, e) 

where G(0, e) = and G(u, 0) = is an analytic function. It is not necessary to assume 
that the higher order terms in G depend on v also, see e.g. jS], but even with such a 
dependence a slightly modified argument would work. As already explained, it is now 
assumed that G(u, e) is such that the invariant curves near the origin for e > are 
compact. 

Hamiltons equation for u reads u = v as before. The period is obtained by solving 
H(u, v;e) = h for v = v(u; h, e) and then by integrating 



T(h,e) = j 



du 



v(u; h, e) 



The idea is to split the integral into two parts; one part near the origin, where the main 
contribution originates, and the rest, which is called T3. In addition the singular integral 
near the origin is split again into two parts, T% which has the same integrand as in the 
previous calculation with G = 0, and a correction T2 which contains the contribution from 
G. The integral T\ will be the most singular, T2 is mildly singular, and T3 is regular. For 
sufficiently small e then T\ dominates and the previous result is recovered. 

To achieve the splitting into T\ and T2 the multiplicative structure of the original 
inegrand has to be recovered. Solving H = h for v and inserting into u gives 

u 2 = v(u; h, e) 2 = Qo(u; h, s)Q(u; h, e) 

where Q = 1 + 0(u) when \u\ < C < 1 for some fixed constant C, and Qq is a polynomial 
of degree 3 in u whose zeroes approach those of the original case with G = when e — > 0. 
Denote by 7c the part of H (u, v ; e) = h for which \u\ < C, and 7c the rest of the invariant 
torus. Then the period integral can be split as 

T(/l ' e) = / im m + / /or h\ (28) 
he \/Q{u,£,h) Jy c ^Q{u,e,h) 

— 1/2 

du + I Q(u h,e) -1^ r du (2g) 



7c \jQa{u,e,h) J lc y/Q (u,e,h) Jj C \/Q{u,e,h) 

=Ti + T 2 + T 3 . (30) 

Now T3 is regular, and gives a bounded contribution, so we can ignore it. The integral T\ 
is singular in the limit h, £ — > 0. T2 is less singular, and in particular bounded, because 
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£ 


1 10 
U£ ' 




0.1 


0.027 


0.28314 


0.01 


0.328 


0.17278 


0.001 


0.478 


0.14656 


0.0001 


0.528 


0.14008 


0.00001 


0.545 


0.13823 


0.000001 


0.550 


0.13767 



Table 1: Numerically measured position u and rotation number of the twistless curve 
for the three times iterated Henon map with parameter k = 1 — £. 

the numerator goes to zero in this limit. So we only need to show that T\ approaches the 
complete integral when e — ► 0. The integral is 



T 1 (h,e) = ^^F(cf ) ,k), ^ = 2arctan^_±^, 

where uq is the real root, ( is the complex roots of Qo, and r and k are as before, 

H + |C| 2 , 2k 2 = l + ^. 

The integral T\ does not diverge because of the modulus k 2 — > 1, but because in the 
prefactor r — > 0. In fact, we already observed that the modulus k 2 < (s/3 + 2)/4 < 1 
inside the first quadrant e > 0, h > 0. The behaviour of r for small e is obtained from 

r 2 = ^u^l + a 2 = D{ 1 )e 1 l A . 

This shows that <p — > it, and the integral approaches the complete integral with the same 
modulus k(j) as before. Therefore even though T\ is an incomplete integral, in the limit 
of small £ its value approaches that of the complete integral T(h, e), which diverges. Since 
the other integrals T2 and T3 are finite the analysis obtained from the complete integral 
over the non-compact curve of H(u, v; e) with G = therefore correctly describes the 
behaviour of the rotation number of the compact invariant curve obtained when G / 0. 

5 Example: Henon Map 

The Henon map in the area preserving case, 
0',y') = (y ~ k + x 2 , -x) , 

illustrates the above. There are saddle-centre bifurcations in the Henon map for which the 
invariant tori for e > are not compact, and hence the vanishing twist cannot be observed. 
This applies to the initial bifurcation at k = — 1 that creates the pair of fixed points, and 
also to many saddle-centre bifurcations that occur for k > 4. However for k = 1 a pair of 
period 3 orbits is created at the origin in a saddle-centre bifurcation of the third iterate 
of the map, and the corresponding invariant tori are compact. One of the period three 
points is located on the symmetry line x = —y. In new coordinates (u, v) = (x, y + x) the 
third iterate of the map expanded near the origin with parameter k = 1 — e is 

( u \ v ') = [(u - v )(l + 4v),v + £ + (u 2 - 2uv + 3v 2 )) + 0(u 3 , v 3 , vlef) • 
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The (exact) location of the fixed points is (u,v) = (±y/—e, 0), with trace of the Jacobian 
2 =F 2-/ z e + 8s =F 8(-e) 3//2 so that the (approximate) multiplier is fi ss 1 + i\/2(-e) 1 / 4 
and /i for the fixed point at u = y/—e. The corresponding rotation number uj of the 
elliptic fixed point is obtained from fi = exp(27riu;), so that u = (— e) 1 / 4 /(v / 27r). For small 
positive e the Henon map possesses compact invariant curves near the origin. A higher 
order normal form would give a G that describes these invariant curves. The Henon map 
is non-integrable, so that only sufficiently irrational invariant curves of the (high order) 
normal form will exist in the Henon map. For the situation under consideration numerical 
experiments show that many of these invariant curves do exist. 

Considering the third iterate of the Henon map turns the pair of period three orbits 
into three pairs of fixed points with heteroclinic connections. In the normal form there is 
only one pair of fixed points, and the unstable fixed point has a homoclinic connection. 
To match the prediction in the case of more than one unstable fixed point the period and 
hence rotation number must be calculated for the heteroclinic connection. For e > there 
are invariant tori on which the dynamics becomes slow near the three points that are 
close to the three bifurcation points of the third iterate of the map. Hence the rotation 
number for the third iterate of the map between two successive such points gives the 
correct rotation number. The results together with the position of the twistless curve are 
shown in table E The values shown converge to the predicted values given in (|15|) and 
(|16j). however, fairly small e are needed to see this. The convergence to the true value 
0.1374244... occurs approximately 
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